8.6 Linearization of Nonlinear Systems 

In this section we show how to perform linearization of systems described by 
nonlinear differential equations. The procedure introduced is based on the Taylor 
series expansion and on knowledge of nominal system trajectories and nominal 
system inputs. 

We will start with a simple scalar first-order nonlinear dynamic system 

x(t) = F(x(t),f(t)), x(t 0 ) given 

Assume that under usual working circumstances this system operates along the 
trajectory x n (t) while it is driven by the system input We call x n (t) and 

respectively, the nominal system trajectory and the nominal system input. 
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On the nominal trajectory the following differential equation is satisfied 

X n (t) = T(x n (t),U n (t)) 

Assume that the motion of the nonlinear system is in the neighborhood of the 
nominal system trajectory, that is 

x(t) = x n (t ) + A x{t) 

where A x(t) represents a small quantity. It is natural to assume that the system 
motion in close proximity to the nominal trajectory will be sustained by a system 
input which is obtained by adding a small quantity to the nominal system input 

f(t) = f n (t) + A f(t) 
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For the system motion in close proximity to the nominal trajectory, we have 


x n {t) + A x(t) = F(x n (t) + A x{t), fn{t) + A f(t)) 


Since A x{t) and A f{t) are small quantities, the right-hand side can be expanded 
into a Taylor series about the nominal system trajectory and input, which produces 


dJ= dJ= 

X n {t) + A X{t) = ^(Xni fn ) + —(x n , f n )Ax(t) + —— (x n , f n )Au(t) 

OX on 


+ higher-order terms 


Canceling higher-order terms (which contain very small quantities 
Ax 2 , A/ 2 , AxAf , Ax 3 ,.. .), the linear differential equation is obtained 


A x(t) = 


dJ= 

dx 


(Xni f n )Ax(t) + 


dJ= 

~df 


(X n , fn)Af(t) 


The slides contain the copyrighted material from Linear Dynamic Systems and Signals, Prentice Hall 2003. Prepared by Professor Zoran Gajic 8-85 



The partial derivatives in the linearization procedure are evaluated at the nominal 


points. Introducing the notation 

dJ= 

Cto(^) = ~ (Xni /n)? 

OX 


bo — o O fn) 

df 


the linearized system can be represented as 


A x(t) + a 0 (t)Ax(t) = b 0 (t)Af(t) 


In general, the obtained linear system is time varying. Since in this course we 
study only time invariant systems, we will consider only those examples for which 
the linearization procedure produces time invariant systems. It remains to find the 
initial condition for the linearized system, which can be obtained from 

Ax(t 0 ) = x(t 0 ) — x n (t 0 ) 
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Similarly, we can linearize the second-order nonlinear dynamic system 
x = t{x, x(t 0 ), x(t 0 ) given 

by assuming that 

x(t) = x n {t) + A x(t), x(t) = x n (t) + A x{t) 
fit) = fn(t) + A/(£), f(t) = /„(*) + A/(J) 

and expanding 

x n + A® = f(x n + A® n , + A®, fn + A/, fn + a/) 

♦ 

into a Taylor series about nominal points x n , x n , / n , f n , which leads to 
A x(t) + a\ Ax(t) + aoA®(f) = b\Af(t) + boAf(t) 
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where the corresponding coefficients are evaluated at the nominal points as 


dT ( . . x dJ= 

ft 1 — 777 ~ ( x ni x n , f n , f n ), CtQ — X, ( x n ? fn 5 fn 


dx 


dx 


dJ= ( . . X dJF/ 

^1 — y x m x n ? /n? fn J ? ^0 — ~df \ ™ f‘ 


The initial conditions for the second-order linearized system are obtained from 


Aa?(f 0 ) = *(fo) — *n(fo), A£(f 0 ) = £(f 0 ) — £ n (f 0 ) 


Example 8.15: The mathematical model of a stick-balancing problem is 

0(f) = sin 0(f) — /(f) cos 0(f) = JF(0(f),/(f)) 

where /(f) is the horizontal force of a finger and 0(f) represents the stick’s angular 
displacement from the vertical. 
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This second-order dynamic system is linearized at the nominal points 


6 n {t) = 6 n (t) = 0, f n (t) = 0, producing 


o\ = 


dT 

x = 0, 


d0 


«o = — 


dT 

~d0 


= — (cos# + fsinO).e„(t)=o = 

1 /»(*) =0 


-1 


, dT 

6l = -r = 0, 

df 


dF\ 

bo = ( -) = —(cos 0)i = —1 


The linearized equation is given by 


m - e(t) = -f(t) 


Note that A 6(t) = 6(t ), A/(£) = /(£) since (9 n (t) = 0, / n (£) = 0. It is 
important to point out that the same linearized model could have been obtained by 
setting sin Q(t) ~ 0(f), cos 0(t) ~ 1, which is valid for small values of 0(t). 
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We can extend the presented linearization procedure to an n-order nonlinear 
dynamic system with one input and one output in a straightforward way. However, 
for multi-input multi-output systems this procedure becomes cumbersome. Using 
the state space model, the linearization procedure for the multi-input multi-output 
case is simplified. 

Consider now the general nonlinear dynamic control system in matrix form 

^x(f) = JF(x(£),f(£)), x(£ 0 ) given 

at 

where x(£), f(£), and T are, respectively, the ro-dimensional system state space 
vector, the r-dimensional input vector, and the ro-dimensional vector function. 
Assume that the nominal (operating) system trajectory x n is known and that 
the nominal system input that keeps the system on the nominal trajectory is given 
by f n (t). 
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Using the same logic as for the scalar case, we can assume that the actual 
system dynamics in the immediate proximity of the system nominal trajectories can 
be approximated by the first terms of the Taylor series. That is, starting with 

x(f) = x n (t) + Ax(f), f (t) = f n(t) + Af (t) 


and 


d 

—Xn(t) = ^ r (x n (t),f n (t)) 

at 


we expand the right-hand side into the Taylor series as follows 


d d 

-Xn + -Ax = ^ + Ax, f n + Af) 
at at 


= f(x n ,f n ) + 


dJ= 

^X J | *n(t) 

t'n(t) 


Ax + 


dJ= 

9f J | x n(0 

t'ri(t-) 


A f + higher-order terms 
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Higher-order terms contain at least quadratic quantities of Ax and Af . Since Ax 
and Af are small their squares are even smaller, and hence the high-order terms 
can be neglected. Neglecting higher-order terms, an approximation is obtained 
d , x fdF\ , x fdF\ 

* Ax(t)= (ftr) | ».(.. Ax(t)+ U 

■f nit) 7 If „{t) 

The partial derivatives represent the Jacobian matrices given by 



( 



__ j^nXn 

x n (t) 
f n(t) 


v9J= x 

dj= x 


9J =i -| 

dxi 

dx 2 

• • • 

• • • n 

C/d&n 

9T 2 



9J=2 

dx\ 

• • • 

• • • 

• • • rj 

OX n 



9J=i 


• • • 

• • • 

9xj 

• • • • • • 

• • • 

9T n 

• • • 

8J= n 

• • • 

• • • • • • 

9T n 

_ dx\ 

dx 2 

• • • 

• • • rj 

9x n _ 


I x n (t) 
fn(*) 
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( 



I x n (t) 
fn (t) 


_ gnXr 


rdT x 

dTi 


dTx -I 

dfi 
dT 2 

df 2 

9 • • 

* * * df r 

dT 2 

dfi 

9 9 9 

• • • 

dT % 

* * * df r 

9 • • 

• • • 

dfi 

9 9 9 999 

• • • 

dT n 

• • • 

dT n 

• • • 

9 9 9 9 9 9 

dTn 

L dfi 

dfi 

• • • 

• • • dfr J 


Note that the Jacobian matrices have to be evaluated at the nominal points, that 
is, at x n (<) and f n (t). With this notation, the linearized system has the form 
d 

— Ax(t) = AAx(f) + BAu(t), Ax(f 0 ) = x(f 0 ) — x„(f 0 ) 

at 

The output of a nonlinear system satisfies a nonlinear algebraic equation, that is 

y(t) = 2(x(t),f(T)) 
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This equation can also be linearized by expanding its right-hand side into a 
Taylor series about nominal points x n (t) and f n (t). This leads to 

y„ + Ay = + (|0 Ax + (^) t) Af 

■f nW !fnW 

+ higher-order terms 

Note that y n cancels term C/(x n ,f n ). By neglecting higher-order terms, the 
linearized part of the output equation is given by 

Ay (t) = CA x(t) + DAf (t) 

where the Jacobian matrices C and D satisfy 
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QpXn _ 


D pXr = 


ag_ 

dx 


| x n (t) 

fn (t) 


rog 1 

50i 


eg 11 

dx\ 

eg 2 

dx 2 

• • • 

• • • o 

9x n 

og 2 

dx\ 

• • • 

• • • 

dQi 

• • • r\ 

OXn 

• • • 

• • • 

dxj 

• • • • • • 

• • • 

dQ p 

• • • 

dQ p 

• • • 

• It • • • 

dg p 

_dx\ 

dx 2 

• • • 

* * * dx n J 


og 

df 


x n (<) 


rdg 1 

501 


50i I 

5/l 

5/2 

• • • 

* * * 9fr 

502 

5/i 

• • • 

• • • 

50i 

* * * 9f r 

• • • 

• • • 


• • • • • • 

• • • 

og p 

• • • 

50 p 

• • • 

• • • • • • 

50 p 

Idh 

5/2 

• • • 

* * * 5/ J 
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Example 8.16: Let a nonlinear system be represented by 

— 77 -^ = xi(t) sinx 2 (t) + x 2 (t)f(t) 
at 

d ^l = Xl (t)e-^+f(t) 

at 

y(t) = 2x\{t)x 2 {t) + x\(t) 


Assume that the values for the system nominal trajectories and input are known and 


given by xi n ,x 2n , and f n . The linearized state space equation of this nonlinear 


system is obtained as 


Ax\ (t) 


sin x 2n xin cos x 2n + f n 


Ax\(t) 

i 

X 2 n 

Ax 2 (t) 


e~ X2n —xi n e~ X2n 


Ax 2 (t) 


2f„ 


a m 


A y{t) = [2x 2n 


2x ln + 2x 2n 


Aaq (t) 
Ax 2 (t) 


Having obtained the solution of this linearized system under the given system input 
A/(<). the corresponding approximation of the nonlinear system trajectories is 
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x n (f) + A x(f) 


Xln(t ) 

_l_ 

Ax\(t) 

X 2 n (t) 


Ax 2 (t) 


Example 8.17: Consider the mathematical model of a single-link robotic ma¬ 
nipulator with a flexible joint given by 

10\{t) + mgl sin0i(f) + k(0\(t) — 0 2 (£)) = 0 

J0 2 (t) - k{p x (t) - 0 2 (t)) = f(t) 

where 0\{t ), 0 2 (t) are angular positions, /, J are moments of inertia, m and / are, 
respectively, link mass and length, and k is the link spring constant. Introducing 
the change of variables as 

xi(t) = 0i (f), x 2 (t) = 0i(£), * 3 (^) = 02 (^), x^(t} = 0 2 (f) 
the manipulator’s state space nonlinear model is given by 
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Xl (t) = X 2 (t) 
jYlcil Jx, 

x 2 (t) =- — sinx 4 (t) — j(xi(t) — x 3 (t)) 


x 3 {t) = x 4 (t) 

X4(t) = j(x! (t) - X 3 (t)) + jfit) 


Take the nominal points as (a?i n , x 2n , x 3n , x 4n , /„), then the matrices A and B 
are 


0 

1 

0 

o" 


"o" 

k+mgl cos xi n 

0 

0 

k 

0 

1 


0 

0 

I 

0 

I 

0 

, B = 

k 

J 

0 

k 

J 

0 


1 

_J _ 


Assuming that the output variable is equal to the link’s angular position, that is 
y(t) = x 4 (t), the matrices C and D are given by 


C = [1 0 0 0], D = 0 
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